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Abstract 

There is a close relation between the problems of second layer nucleation in 
epitaxial crystal growth and chemical surface reactions, such as hydrogen re- 
combination, on interstellar dust grains. In both cases standard rate equation 
analysis has been found to fail because the process takes place in a confined 
geometry. Using scaling arguments developed in the context of second layer 
nucleation, I present a simple derivation of the hydrogen recombination rate 
for small and large grains. I clarify the reasons for the failure of rate equations 
for small grains, and point out a logarithmic correction to the reaction rate 
when the reaction is limited by the desorption of hydrogen atoms (the second 

order reaction regime). 
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Rate equations are a widely used tool for the mathematical description of kinetic phe- 
nomena. Despite successful applications in many different areas of science, however, there 
are instances where their predictions are known to be qualitatively incorrect. A prominent 
example are diffusion-controlled reactions in low dimensions [1]. The purpose of this Rapid 
Communication is to point out the close relation between two kinetic problems, second layer 
nucleation in epitaxial crystal growth [2-5] and chemical reactions on dust grain surfaces 
in interstellar space [6-9], in which deviations from rate equation theory arise because the 
processes of interest take place in a spatially confined geometry. 

In both cases one is concerned with the rate of a reaction confined to a finite surface 
region, onto which particles impinge and from which they leave again if no reaction takes 
place. In the case of second layer nucleation the confining region is the top terrace of a 
two-dimensional island, bounded by a step edge, while in the astrophysical context it is the 
(closed) surface of a small grain. In the first case the atoms leave the island by crossing the 
step, while in the second case they desorb from the grain surface^. 

The two problems also differ in that in second layer nucleation only the first nucleation 
event on the island is relevant - all atoms deposited after this event merely contribute to 
the growth of the second layer nucleus - whereas in the problem of grain surface reactions 
the reactants are usually assumed to desorb from the grain once they have formed, so that 
a steady state with a constant production rate is established. In this sense, second layer 
nucleation is a first passage problem [5], while grain surface kinetics deals with stationary 
processes. 

For both problems it has been found that standard rate equation analysis fails when 
the mean number of particles present on the surface becomes small compared to unity (the 
lonely adatom limit of Ref. [3]). In this fluctuation-dominated regime the reaction becomes 



^At high temperatures desorption plays a role also in epitaxial growth; see [10] for an analysis of 
two-dimensional nucleation in the presence of desorption. 
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a rare event which has to be treated statistically. Here I show how the time scale analysis 
developed in the context of second layer nucleation [3,4] can be applied to the problem of 
hydrogen recombination, one of the most important grain surface reactions of astrophysical 
interest. In this reaction two hydrogen atoms meet to form a H2 molecule, in analogy to 
the formation of a second layer nucleus in the case where dimers are stable (irreversible 
nucleation, critical nucleus size i* = 1). I provide a simple derivation of the recombination 
rate in all three reaction regimes. I then discuss the reasons for the breakdown of the rate 
equation theory for small grains, and I point out the existence of a logarithmic correction 
to the recombination rate in the second order reaction regime, which was overlooked so far 
in the astrophysical context. 

The notation follows that of [8]. Atoms impinge on the grain surface at rate / per 
adsorption site, they hop at rate a and desorb at rate W. The number of adsorption sites 
on the grain is S, and the quantities F — fS and A — a/S refer to the whole grain. To 
make contact with the notation of [3,4] we introduce the lifetime r of an atom on the grain, 

T = 1/W, (1) 

the diffusion time td required for an atom to explore the whole grain, 

td^ S/a=l/A, (2) 

and the mean time interval At between the arrival of two atoms on the grain. 

At = 1/F. (3) 

The estimate (2) ignores a logarithmic correction arising from the distinction between the 
total number of sites visited by an adatom during its lifetime, A/'aii(r), and the number 
•^dis{T) of distinct sites visited; we will return to this point towards the end of the paper. 
Following [8] we introduce the dimensionless quantities 

Svisit = a/W = Sir/ro) (4) 

and 
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s^,^,^w/f^s{/:^t/T). (5) 

For large grains {S > Syisit and 5" > Svacant, respectively) Syisit is the number of sites visited 
by an atom during its lifetime, while Syacant is the number of vacant sites surrounding a site 
occupied by an atom. If r > td the atom has time to eplore the whole grain, so that the 
number of sites visited becomes equal to S, while for t < At the grain is typically occupied 
by zero or one atoms, so that an atom is surrounded by S vacant sites; note that the mean 
number of atoms on the grain is equal to 

n = r/At = F/W. (6) 

The relative ordering of the three dimensionless numbers S, s visit and s vacant defines six 
different regimes. It turns out, however, that the behavior of the reaction rate in a given 
regime depends only on which of the three numbers is the smallest. Thus we are left with 
three different regimes, which will be discussed in the following. 

The fluctuation-dominated (lonely adatom) regime is characterized by the condition [8] 

S -C min [Svisit , ^vacant] ■ (7) 

According to (4,5) this is equivalent to (i) r <^ At and (ii) r ^ r^. The first condition 
simply implies that the mean number of atoms on the grain (6) is small compared to unity. 
In this limit, the probability that a freshly deposited atom participates in a recombination 
event is equal to the probability Pi that a second atom is already present on the grain, times 
the encounter probability Penc that the two atoms, once both present on the grain, will meet 
before one of them redesorbs. The first probability is given by [3] Pi — r/ (r + At) ^ r/At — 
F/W . Multiplying this with the number of atoms F = 1/ At arriving on the grain per unit 
time, we find that the total recombination rate (the number of molecules produced by the 
grain per unit time) is given by 

whenever n <^ 1. The second condition t ^ td implied by (7) means that, if two atoms 
are present simultaneously on the grain, they will meet with unit probability, because each 
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atom spends enough time on the grain to visit all sites many times. Thus we have Penc = 1 
and 

R = i^lonely = T / [Mf = /W. (9) 

In the small grain limit where (7) is satisfied, this expression is exact. The recombination 
efiiciency defined by 77 = R/{F/2) is then equal to 2F/H^, in agreement with the result 
obtained in [8] from the full solution of the master equation. In second layer nucleation, the 
lonely adatom regime is realized under typical experimental conditions when i* = 1 [3]. 

In the remaining two regimes the finite grain size is irrelevant. These regimes can there- 
fore be derived from standard rate equations [11]. According to the scaling of the reaction 
rate with hydrogen flux, these two regimes will be referred to as the first and second order 
regimes (note that the lonely adatom regime is also of second order in this sense, see (9)). 

The first order regime is characterized by the condition 

Svacant < mm[S, Syisit] . ( 10) 

In this regime the mean number of atoms on the grain is large compared to unity (r ^ At) , 
and the number of sites visited by an atom during its lifetime (= minfjS", Syisit]) is large 
compared to the number of vacant sites surrounding each atom. Hence each freshly deposited 
atom will certainly find a partner. In this case the recombination rate is 

R^R^^ F/2 (11) 

and the recombination efficiency is equal to unity. In the context of second layer nucleation, 
this regime corresponds to the nonstationary scaling regime I of [2,4], where the energy 
barriers at the step edges are so strong that r S> At, and nucleation occurs as soon as a 
second atom is deposited onto the island. 

It is easy to derive an interpolation formula connecting (9) and (11) under the assumption 
that the diffusion of the reactants, as well as the dcsorption of the hydrogen molecule, is 
very fast. Then the only possible values for the number of atoms on the grain is or 1, 
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since recombination occurs instantaneously as soon as a second atom is deposited. The 
probability Pi for having one atom on the grain satisfies the master equation 

dPi/dt = FPo - (F + W)Pi = F(l - Pi) - (F + W)Pi (12) 

with the stationary solution Pi — F/{W + 2F). Since Penc — 1, the recombination rate is 

The remaining, second order regime is characterized by the condition 

Svisit < min [S, Syacant] • ( 14) 

The condition Syisit ^ S (i.e., r ^ To) has no clear analogue in the context of second layer 
nucleation, because there the lifetime of the atom is always at least as large as the diffusion 
time (the atom has to reach the island edge before it can leave the island). It is useful to 
distinguish the cases S ^ Svacant (large grains) and S <C Svacant (small grains). In the first 
case many atoms are present on the grain simultaneously. The encounter probability for a 
freshly deposited atom to find a partner during its lifetime is equal to Penc = Svisit/svacant- 
Thus the recombination rate can be estimated as 

fi = fl, = fp.„ = F_ = ^ (15) 

and the recombination efficiency is rj — 2af/W'^. For small grains the probabihty that a 
freshly deposited atom will undergo a reaction is the product of the probability Pi = F/W 
that a second atom is present on the grain when it arrives, and the encounter probability, 
which now reads Penc — Syisn/S (note that the second atom occupies any site of the grain 
with equal probability). The recombination rate then becomes 

^=Ty^'enc = -^, (16) 

which is identical to (15). Thus the second order result hold irrespective of whether the 
number of atoms on the grain is large compared to unity or not. This is similar to certain 
regimes of second layer nucleation with critical nucleus size i* > 1 [4] . 
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The last observation shows that the failure of rate equation theory cannot generally be 
attributed to the fact that the number of reactants on the grain is small compared to unity. 
To clarify this issue, it is instructive to consider the amount by which the reaction rate is 
overestimated by the rate equation expressions when the condition (7) is satisfied. Two 
cases have to be distinguished. For Syacant ^ Syisu we have to compare (9) to Ri, which 
yields 

i?lonely/i?l = F/W = r/At. (17) 

This is precisely the mean number of atoms n on the grain. Thus in this case it is justified 

to say that the rate equations fail because n -C 1. On the other hand, for Syisit ^ -Svacant we 
have 

i?lonely/i?2 = SW/a = ^/sw (18) 

This ratio is not related to the mean number of atoms on the grain. Instead, it is equal 
to the ratio between the number of distinct sites visited by an adatom during its lifetime 
(■A/dis — S) and the number of all sites visited (jVaii — a/VF = Syisit)- Castellano and Pohti 
have emphasized that it is the distinction between A/dis and A/lu which is responsible for the 
failure of rate equation theory in the case of second layer nucleation [5] . 

In the second order regime this distinction implies a logarithmic correction to the expres- 
sion (15). The number of distinct sites visited by a two-dimensional random walk during 
time r is given by Adis('?") ~ 7rar/ln(ar) [12]. The diffusion time should therefore be de- 
termined more accurately by AfdisiTo) — S, and the expression (4) has to be replaced by 
Svisit — A/'dis(T)- Using this in the derivation of the reaction rate in the second order regime 
we find 

J. _PS_ npSa 

The experimental estimates for the activation energies for H diffusion and desorption 
on amorphous carbon surfaces quoted in [8] show that the logarithmic factor in (19) is 
ln(a/H^) w 15 at a temperature of 10 K. 



To go beyond the order of magnitude estimate (19) would require a detailed analysis of 
the random walk problem on the grain surface along the lines of [5]. The master equation 
analysis of [6-8] is not sufficient for this purpose, because it keeps track only of the total 
number of atoms on the grain, thus ignoring the spatial aspects of the encounters between 
the atoms. As a first step in this direction, I present here a simple model calculation of the 
encounter probability Penc for the case of small grains (in the sense of n <C 1), which allows 
to interpolate between the expressions (9) and (19); a similar approach has previously been 
applied in the context of second layer nucleation [13]. 

We take the grain surface to be a circular disk of radius L. One of the reacting atoms 
is assumed to be immobile, and is placed at the center of the large disk in the form of a 
small disk of radius I. In this way the problem becomes radially symmetric. The deposition, 
diffusion and desorption of the second atom is described by the stationary diffusion equation 

where n(r) is the probability density for the second atom in polar coordinates, / is the 
impingement flux per unit area, and the diffusion coefficient is D = aP. The boundary 
conditions for (20) are reflecting at the outer boundary, dn/dr — at r — L, and absorbing 
at the inner boundary {n{r — I) — 0). The encounter probability is then obtained as the 
ratio of the diffusion flux incident onto the inner boundary to the total deposition flux 
incident onto the surface. 



2nlD dn/dr\ 



The solution to this diffusion problem can be found in the thin film growth literature [14] . In 
the limit where the atom radius / becomes the smallest length in the problem, the encounter 
probability takes the form 

Penc = ^ (ln(a/W) + 2K,{^SW/a)/h{^SW/a)y' , (22) 

where /i and Ki are incomplete Bessel functions. For small arguments Ki[x) / Ii{x) 

and Penc — 1, while for large arguments Ki{x)/Ii{x) and (22) reduces to Penc ~ 
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4:a/{SWla{a/W)). Inserting the second result into (8) we obtain an expression identical to 
(19), apart from the factor tt which is replaced by a factor 4. 

To summarize, I have shown in this paper how recent results from the theory of second 
layer nucleation can be applied to the closely related problem of hydrogen recombination on 
dust grain surfaces. The expressions (9), (11) and (15) for the recombination rate have been 
obtained before [8], but the present derivation seems more straightforward. The analysis 
shows that the failure of the rate equation approach is not always related to the fact that the 
mean number of atoms n on the grain is small compared to unity. The logarithmic factor 
in the second order regime, Eq.(19), and the interpolation formulae (13) and (22) are new. 
The derivation of the encounter probability for realistic grain geometries using the methods 
of [5] seems like a promising problem for future studies. 
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